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Abstract 

We study a one-dimensional model of free fermions with gl(l|l) supersymmetry and demonstrate 
how non-diagonal boundary conditions can be incorporated into the framework of the graded Quan- 
tum Inverse Scattering Method (gQISM) by means of super matrices with entries from a superal- 
gebra. For super hermitian twists and open boundary conditions subject to a certain constraint, 
we solve the eigenvalue problem for the super transfermatrix by means of the graded algebraic 
Bethe ansatz technique (gABA) starting from a fermionic coherent state. For generic boundary 
conditions the algebraic Bethe ansatz can not be applied. In this case the spectrum of the super 
transfer matrix is obtained from a functional relation. 
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I. INTRODUCTION 



For a long time studies of quantum integrable models in one spatial dimension have led 
to important insights into the properties of many body systems and provided a sound basis 
for the understanding of the non perturbative phenomena which arise due to the interplay 
of interactions and strong quantum fluctuations in low dimensional systems (see e.g. A 
special way to introduce free parameters into these systems is by variation of their boundary 
conditions. Considering all possible classes compatible with the integrability allows for a 
complete classification of their low-energy quantum critical behaviour on one hand but 
also to study in detail the effect of embedded impurities and contacts to an environment. 
Recently, there has been increased interest in twisted or non-diagonal boundary conditions 
which break certain bulk symmetries of integrable quantum spin chains 
although their hamiltonian is a member of a commuting family of operators the established 
algebraic schemes for the computation of the spectrum fail unless additional constraints to 
the boundary conditions are in place. For spin 1/2 chains there has been some progress 
using functional methods, but quite a few open questions remain. Even less is known for 

though integrable non-diagonal 



m etry . r 



121 ] the solution of the spectral 



open boundary conditions have been constructed [9|, LLCJ, 111 
problem is restricted to diagonal ones so far. 

In this paper we study this problem for the simplest possible case of spin chains with 
$j[(l|l) supersymmetry. Since the corresponding bulk system describes free spinless fermions 
on a lattice this should provide a toy model to investigate in particular the applicability of 
functional methods to the solution of the spectral problem. We begin with a short review 
of the graded Quantum Inverse Scattering Method [l3,Q, [l^. Using a Grassmann valued 
super matrix representation of the Yang Baxter algebra, spin chains subject to twisted pe- 
riodic boundary conditions can be embedded into this framework and are solved exactly. In 
Section HTT1 we construct the jj[(1|1) su per spin chain with generic open boundary conditions 
based on Sklyanin's reflection algebra [16] . We study the spectrum of these super spin chains 
for certain classes of reflection matrices using the algebraic Bethe ansatz and finally extend 
this solution to generic boundaries using functional methods. 
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II. GRADED QUANTUM INVERSE SCATTERING METHOD 



The fundamental objects considered within the framework of the graded Quantum Inverse 
Scattering Method (gQISM) are representations T(v) of the graded Yang-Baxter algebra 
(gYBA) 

Ru(u — v)T («) T iy)=T (v) T (u)R 12 (u - v) . (2.1) 

The indices 1 and 2 label the linear spaces Vi j2 into which the respective operators are 
embedded by means of the super tensor product <8> s , defined through 

(A ® s B)(C ® s D) = (-lyPM^AC ® s BD , (2.2) 

where p(X) refers to the parity function defined in the appendix. That is, to be precise 

T (u) = T{u) ® s l, T{u) = t® s T{u) , 

(2.3) 

R 12 (u) = R(u) ® s 1 , R 23 = l® s R(u) and R n (u) = P 23 R 12 (u) P 23 ■ 

Here Pij is the graded permutation operator that interchanges two spaces Vi and Vj according 
to P(x ® s y) = (— iy( x )P(y>(y (g) s x). The i?-matrix is subject to the consistency condition 

R 12 (u - v)R 13 (u)R 23 (v) = R 23 (v)R 13 (u)R 12 (u - v) , (2.4) 

known as Yang-Baxter equation (YBE). As a consequence one obtains local representations 
Loj(u) = Roj(u) of the gYBA by a graded embedding of the i?-matrix. These Lax-operators 
Loj(u) act on an auxiliary space Vo, whereas their entries act on the j-th quantum space 
Vj. Due to its comultiplication property, the gYBA allows for the construction of global 
representations as products of Lax-operators. This results in a particular representation on 
the auxiliary space and the tensor product of the quantum spaces V q — V\ ® s V 2 ® s • • • ® s Vn, 
the monodromy matrix 

T(u) = L m (u)L Q)N _i(u) . . . L 01 (u) . (2.5) 
Taking the supertrace flAlip of this monodromy matrix, we obtain the super trans fermatrix 



t(u) = str{T(w)} which generates a set of commuting operators on V q . In particular, 
it is related to an integrable hamiltonian with periodic boundary conditions defined by 
H = d u \nT(u)\ u=0 . 

For the g[(l|l) supersymmetric representations of the gYBA considered here, this con- 
struction leads to a model of free spinless fermions on a one-dimensional lattice with N sites. 



15|]) is 



In the case of periodic boundary conditions the hamiltonian reads 

N 
3=1 

The corresponding i?-matrix (cf. 

fu + 1 \ 
u 1 
1 u 

\ «-v 

and a graded embedding yields 



R{u) 



rv R(u) = PR{u) 



fl + u \ 
1 u 
It 1 

v i-v 



L 0j (m) = -Roj(m) 



u + ej :1 ej,2 



U + Uj Cj 



Cj U — 7lj 



(2.6) 



(2.7) 



(2.8) 



Generally we will define the hamiltonian density in terms of the checked i?-matrix via Hij 

d u Rij{u)\ u= Q. 

A. Super hermitian twists 



The simplest generalization of periodic boundary conditions are twists. They can easily 
be incorporated into the above scheme by making use of the comultiplication property again. 
Let the twist matrix K be a representation of the gYBA on the auxiliary space. Then K-T{u) 
is another global representation producing the super transfermatrix 



t(u) = str { K ■ T(u) } = str { K L 0A r(M)L ,iv_i(M) . . . L 01 (u) } , 
which results in a modified hamiltonian on V q which contains a boundary term 

N-l 

As a specific twist matrix we choose 



a dS . 

K= I a,beR , deC, 

(dSf b 



(2.9) 



(2.10) 



(2.11) 



where £ is the sole generator of CGi (see Appendix IA 2ft . Notice that this is the most 
general CGi super matrix being hermitian with respect to the operation (IA13|) . Taking into 



account the properties of Grassmann numbers, K can be diagonalized by a super unitary 
transformation 

1 (Ua-b) d£* \ 

U= \ rx U ] U = UU ] = 1, (2.12) 

a ~ b y d*£ i(a -b)j 

such that 

K = U ] KU= \ a I . (2.13) 

V ») 

The Lax-operators (12. 8p are super matrices over the algebra described in Appendix IA 11 
hence the comultiplication (12.91) will lead to products between fermionic operators (1A3j) and 
Grassmann numbers. For homogeneous elements C G T and G G CGV we define 

[C, G} ± = and p{CG) = p(GC) = p{G) + p{C) mod 2 . (2.14) 

In the periodic case ( 12.61) . the spectrum can be obtained by means of the graded algebraic 
Bethe ansatz (gABA) with the Fock-vacuum as a reference state. For diagonal (or upper 
triangular) twist matrix K the Fock-vacuum would still provide a suitable reference state 
for the gABA. For more general twists a different pseudo vacuum has to be used. 
Using the cyclicity of the supertrace we rewrite the super transfermatrix (12.91) as 

t(u) = Str | ^Z 0A r(M)L(yV_i(u) . . .L m (u) | , (2.15) 

with transformed Lax-operators 



L 0j {u) = tf L 0j (u)U 



J a—b 3 a—b J ] a—b 

\ ca - ^t-£ b u-rij- — h S^c\ + -^- h S Cj 

\ J a—b J a—b J a—b J / 



(2.16) 

By means of a super unitary transformation on the quantum space Vj the Lax-operator 
Loj(u) can be written in the form (12.81) : setting 

P^—^ rx p* = —,S (2.17) 

a — o a — b 

we define unitary operators 

Qi = l + fa + pc) = e pc ^ Q] = l- pc) - p»c, = e -("4+^) (2. 18) 
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that map the fermionic creation and annihilation operators according to 



Cj = QjCjQj = Cj - p and cj = Q)c)Qj — cj — 



rv rij = cjcj = rij + pcj — p**Cj , = CjCj = 1 — rij = 1 — rij — pCj + p^Cj . 
In terms of these new fermionic creation and annihilation operators we obtain 



L 0j (u) 



u + rij cj 



Cj u — rij 



After this transformation the gABA can be applied with the new Fock vacuum 

1 10> 



1 6) = e -"£f=i4i 



as the reference state. Note that the local Fock vacua 



(2.19) 



(2.20) 



(2.21) 



|6;> = Q]|0j> = |0;>-p|1j> (2.22) 
are fermionic coherent states, i.e. eigenstates of the annihilation operator Cj\0j) = p\0j). 



III. GRADED REFLECTION ALGEBRA 



We will now extend Sk 
open boundary conditions 



yanin's formalism for the treatment of integrable systems with 



16| in a way that makes it applicable to supersymmetric models. 



Following 17|, [l8|, for a given i?-matrix we introduce two associative superalgebras 71 and 
T + , subject to the graded reflection equation 



Rn(u - v) T {u)R 21 {u + v) T (v) 

= T_ (v)R 12 (u + v) 71 (u)R 2 i(u-v) 



and to the dual graded reflection equation 



Rf^( v - u ) T+ stl (u)R 12 (-u - v) T 4 



1st 2 / 



T™ t2 (v)R 21 (-u - v) T + stl (u)Rf 2 list2 (v - u) 



respectively, whereas the new matrices R and R are related to the i?-matrix via 



Rt*(-u-v)R%\Hu + v) 



1-12 V " U J AL 21 

? ist l ( n, n<\ E> ist 2 / 



R^(-u-v)R^(u + v) 



1 and 
1 . 



(3.1) 



(3.2) 



(3.3) 
(3.4) 



Moreover the i?-matrix (12.71) satisfies the unitarity condition R^iu — v)R 2 i(v — u) ~ 1. 
Under these conditions it is possible to show that the super transfermatrices 



T(U 



str { T + (u) T (u) } (3.5) 



provide a family of commuting operators, i.e. [t(u),t(v)] = 0, Wu,v G C. 

Now open boundary conditions can be described by two auxiliary space matrices K_{u) 
and K + {u) satisfying the reflection equations (13.11) and (13. 2L Up to normalization, the 
restriction to CG\ essentially 1 yields solutions 

(a ± b± S \ 
(3-6) 
f±£ 9 -a±J 

with complex coefficients a±,b± and f±. 

Let T{u) be a representation of the gYBA (12.11) . Then T(u)K_(u)T^ 1 (—u) is a further 
representation of the graded reflection algebra 71 and we have 

t(u) = str { K+i^T^K^T-^-u) } . (3.7) 

The i?-matrix is regular, that is -R(O) = P, and for convenience let us choose the normaliza- 
tion such that K_(0) = t. Since K + (0) has a vanishing supertrace we compute the second 
derivative of the super transfermatix (13.71) and - bearing in mind that the i?-matrix (12. 7h 
complies with the unitarity condition only up to normalization - find 

d 2 



du 

with the open chain hamiltionian 



2 r(u) 



= 8[l + a + ]H (3.8) 

u=0 



3=1 



1 d n 



(3.9) 



u=0 

Now we may address the question of what type of boundary terms the matrices K_ and K + 
do generate, i.e. in what way such boundary conditions affect the hamiltonian of the given 
model. Using the expressions (13. 6p explicitly, the hamiltonian (13.91) can be written as 

1 / a_ gL£ \ 1 I a+ d + £ \ 

^s^ +i+ H/-»4 + ^U-J/ (310) 



1 Constant matrices of the form K± = (K±(u) — l)/u can be employed as well. 
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In using standard representations of (1A3j) and by exploiting (1A14j) we can express the two 
matrices from the latter equation by elements of the combined superalgebra. The first matrix 
yields 




S\ 

0/ \ £ 9 





-1 



a_(l - 2n) + d_E c - f.&S 



(3.11) 



(3.12) 



(3.13) 



and after repeating this procedure for the second matrix, the entire hamiltonian reads 

AT-l 

2 

1 



+ 



2a_n 1 + d_£ Cl - fS^c\ 

a + — 2a + riN + d + £cN — f+£^c^ N 



(3.14) 



2(1 + a+) 

We point out that the non-diagonal boundary terms, which do not preserve the particle 
number, are Grassmann valued (i.e. ~ £). Such terms may arise, e.g., in the description 
of the system coupled to a fermionic environment after integrating out the bath degrees of 
freedom. 



IV. GRADED ALGEBRAIC BETHE ANSATZ 

In this section we show how the spectral problem for the hamiltonian f)3.14p can be 
solved by means of a graded algebraic Bethe ansatz. For notational convenience we set 
T{u) = T(M)i^_(M)T _1 (— u) and consider T(u) as a 2 x 2-matrix 

/ (a{u) B(u)\ t 
T(u) = (4.1) 

\C(u) V(u)J 

on the auxiliary space. The reflection equation (13. lft gives commutation relations between 
the quantum space operators A(u), B(u),C(u) and T>(u) of which the following three are of 
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particular interest 
B(u)B(v) l ~ U + V 

A{u)B{v) 



B{v)B{u) 



V{u)B{v) 



1 + u — v 

(1 — u + v)(v + u) 
(1 + u + v)(v — u) 

(1 — u + v){v + u) 



B{v)A{u) + 
B{v)V(u) 



1 + U + V 

1 



B(u) 



U + V 
u — V 



A(v)-V(v) 



U + V 



V{v) - A{v) 



(l + u + v)(v — u) v / v ' ' 1 + u + v lu — V 
Let |0) be a pseudo- vacuum upon which T(u) acts as an upper triangular matrix, i.e. 



(4.2a) 
(4.2b) 
(4.2c) 



T(u)|0) 



(4.3) 



A(u)\0) B(u)\0)\ _ (a(u)\0) * ^ 
C(u)|0> V(u)\0)) ~ \ 6(u)\0) / 
Here a(u) and 5(-u) are scalar functions, called parameters, that are to be determined later 
on. They are eigenvalues to A(u) and V(u) for the eigenstate |0). 



A. Diagonal boundary conditions 



We begin by considering diagonal boundary matrices K_ and K + , i.e. 



Kju) 



1 + wa. 



and K+iu] 



1 + wa_ 



1 — ua 4 



(4.4) 



1 — wa_ 

This yields the super transfermatrix 

t(u) = str { K + {u)T (u) } = (1 + ua+).A(u) - (1 - ua + )V{u) . (4.5) 



Using the commutation relations (14. 2 aft to (I4.2c[) we find B{v\) . . . B(vm)\0) to be an eigen- 
state of t(u) with eigenvalue 

(l — u+ v e )(v e + u) 



A(u) 



n 



L (1 + W + f<)(-U£ - It) 



(1 + ua + )ac(u) — (1 — ua + )5(u) , (4.6) 



provided that the Bethe ansatz equations 

a(vj) 1 — a + Vj 



(4.7) 



5(vj) l + a + vj 

are satisfied. Here the functions a(u) and 5{u) are obtained from the action of T(u) on the 
Fock vacuum |0) 

a{u) B(u) 



T(u)\0)=T(u) K„ (u)T- L (-u)\0) 



5(u) 



|0>. 



(4.8) 



Using (J23D and O we find 



a(u) 
5{u) 



1 — v? 



1 — V? 



N 



(1 + ua^)[u + 1] 



2iV 



1 - ua)u 2N + (1 + m. 



u 



2JV 



1 + 2u 



u + 1 



1 1 



2N 



(4.9) 



- 1 



Therefore the Bethe ansatz equations 



Vj + 1 



N 



1 — a + Vj 1 — aAvj + 1} 



1 + a + {vj + 1) 1 + a-Vj 



(4.10) 



determine the quantization of single particle momenta of the free fermions due to the bound- 
ary conditions. 

Finally, we find an explicit expression for the operators B(u), that generate eigenstates 
of the super transfermatrix: 

j'-i 



B{u) 



, N o N 

u \ 2 



1-u 2u + 



+ 



u 



u + 1 



II 



u + 1 



3-1 



(4.11) 



ci . 



B. Quasi-diagonal boundary conditions 

Application of the graded Bethe ansatz for non-diagonal boundary matrices is only pos- 
sible when a suitable reference state can be found. Here we consider a super hermitian left 
boundary matrix K + 



t\,Jji) I + a | a+ d+S I with a + eR und d+ E C 
d*E^ -a, 



(4.12) 



which is diagonalized by the super unitary transformation 



U 



1 / 2ia+ d + £ s 



2a + \dl£ 2xa A 



K + (u) = U j K + (u)U 



1 + ua A 



1 — ua^ 



(4.13) 



Now we proceed as in Section III Al the transformation U leaves the Lax-operators shape- 
invariant, and we find 



u + fx ■ c • 
L oj (u) = rfL 0j (u)U- 1 ' ' 

Cj u — rij 



(4.14) 
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where Cj = Cj — p and cj = cj — p"; but now we have 



p^^fl ^ p* = ^£. (4.15) 



2a + 2a^ 



Due to the cyclicity of the supertrace the super transfermatrix can be written as 

r{u) = str { K + (u)T(v) } = str | K + {u)f{u) J , (4.16) 

where 

~ , ~ n /.A(u) 

Tu =Tu ^ ur 1 -^ = ~ ~ • (4-17) 

\C(u) V(u)J 

Here we have introduced T(u) = Lqn{u) . . . Lqi(u) and K_ is the transformed right boundary 
matrix 



a. 



1 / a + (l + wa_) (a + d- — d + a)u£ . 
K^u) = U ] K^u)U = — [ +l ; l+ ]• (4-18) 

(a + f- — d* + aJ)uE^ a + (l — wa_) 



Now, choosing the parameters in (14. 18f) to satisfy the constraint 

a+ /_ = d; a _, (4.19) 

the transformed boundary matrix EL is upper triangular and the graded algebraic Bethe 
ansatz can be performed again with a pseudo vacuum constructed from the fermionic coher- 
ent state (I2.2ip by using the definition (14. 15ft for p (see Ref. 2] for a similar approach in the un- 
graded case). Furthermore, since the transformed quantum space operators A(u), B{u),C{u) 
and T>{u) obey the same fundamental commutation relations (I4.2al) to f)4.2cp as their original 
counterparts, the Bethe ansatz equations f)4.10p remain unchanged. 

Compared to the diagonal case we find that the addition of non-diagonal boundary pa- 
rameters subject to the constraint (14. 19f) does not affect the eigenvalues of the super transfer- 
matrix: the energy spectrum of the chain is determined by the diagonal parameters a± of 
the boundary matrices alone. The Bethe states are generated by the action of the operator 
B on the new pseudo vacuum. Due to the unitary transformation it contains a Grassmann 
valued shift 

Nr. N 



B(u) 




(4.20) 

— a_ ut . 
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Therefore, the Bethe states £>(t>i) . . . B(vm)\0) are linear combinations of states with up to 
M particles added to the coherent state Fock vacuum (12.211) . 



C. Generic boundary conditions: functional relations 

Finally, we want to address the question to what extent the spectral problem of the 
g[(l|l)-model can be solved if we choose more general boundary matrices than those allowed 
by the constraint (14.19p . In this case a reference state suitable for the application of the 
gABA is not available. 

In the case of spin 1/2 chains without grading this question has been addressed by ex- 
ploiting certain functional relations obeyed by the eigenvalues of the transfer matrix as a 
consequence of integrability of the model (see e.g. [si, Q, O, Q] ) ■ To obtain such a functional 
relation for the model considered here we begin with the representation ( 14. 6ft of the eigenval- 
ues in terms of roots of the Bethe equations. Note that only the eigenvalues of the boundary 
matrices enter this expression in the cases studied above. 

Let k± 2 be the eigenvalues of the boundary matrices K±(u), then (14. 6p can be rewritten 
as a functional relation for an unknown function q(u) 



By construction A(u) is a polynomial in u. Therefore Eq. (I4.2ip has to be complemented 
with the condition that its RHS is analytic. In particular the residues at the zeroes of the 
unknown function q(u) have to vanish. With a polynomial ansatz 




(4.21) 



where f(u) is a known function: 



f(u) = k + a(u) — k+8(u) 




(4.22) 



M 




(4.23) 



this leads immediately to the Bethe equations (14.101) . 
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For spin 1/2 chains it has been observed [5j, [7|, that the functional equations such as 
f)4.2ip hold both in the case of diagonal or quasi-diagonal and in the generic off-diagonal 
boundary conditions: there, only the eigenvalues of the boundary matrices enter the equation 
explicitly while the deviation from constraints such as (14. 19ft in the non-diagonal case changes 
the asymptotic behaviour of its solution. This leads to non-polynomial solutions q(u) to the 
corresponding difference equations and therefore Bethe like equations are not easily obtained. 

Based on this observation we propose that the eigenvalues of the super transfermatrix 
( 13. 7p satisfy Eq. (I4.2ip with f{u) parametrized by the eigenvalues of the generic boundary 
matrices K±(u) as in (I4.22p . We have verified this hypothesis for small system sizes where 
we are able to explicitly construct the super transfermatrix as a square even super matrix of 
corresponding finite dimension. Taking into account the peculiarities arising from grading 
as well as the nilpotency of Grassmann generators, it is perfectly possible to perform an 
exact diagonalization by the use of computer algebra systems. For chains with up to iV = 6 
sites we have computed the eigenvalues for the most general boundary matrices K^(u) 
and K + (u) and found that the functional equation (I4.2ip is indeed satisfied. Unlike the 
situation for spin 1/2 chains, however, the functions q(u) are still polynomial as in (I4.23P 
which allows to compute the eigenvalues by solving the Bethe equations (14.101) for generic 
boundary conditions! 

As an simple example we consider a system with just one site, i.e. N — 1: the exact 
diagonalization of the corresponding super transfermatrix yields the two eigenvalues 



On the other hand, assuming that the eigenvalues satisfy (I4.2ip with polynomial q(u) (I4.23P 
we can determine the values of the parameters ve from the requirement, that A (it) has 
vanishing residues at the poles at u — v £ and u — — 1 — V£. For M = we immediately 
obtain A + (u) while for M = 1 we find 




2u 



(1 + a + + u{u ± 1) [a + + a_(l + a+)]) . 



(4.24) 




(4.25) 



and thereby recover the second eigenvalue A (u). 
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V. SUMMARY AND CONCLUSION 



In this paper we have studied g[(l|l)-symmetric super chains of free fermions subject to 
generic non-diagonal - in general Grassmann valued - boundary fields breaking the U(l) 
particle number conservation of the bulk system. The boundary conditions could be em- 
bedded into the reflection algebra formalism resulting in quantum integrable models. For 
the solution of the spectral problem we have applied the graded algebraic Bethe ansatz 
for a class of boundary conditions satisfying a constraint (14.191) . In these cases both the 
eigenvalues and the eigenstates of the super transfermatrix are obtained by the action of 
creation operators on a suitably chosen reference state. For generic boundary conditions 
such a vacuum state could not be constructed. Motivated by recent findings for spin chains 
without grading we have proposed the hypothesis that the eigenvalues can still be obtained 
from Bethe equations and verified this conjecture for small system sizes using numerical 
methods. In this case, however, it is not clear how the eigenstates are parametrized by the 
Bethe roots. 

Although the case of g[(l|l)-symmetric super chains is particular simple since the resulting 
hamiltonian describes free particles, our results indicate that it may be easier to deal with 
non-diagonal boundary fields in integrable super spin chains than in models without grading. 
A straight forward extension is to the g-deformation of the system presented here. Non- 
diagonal solutions to the reflection equations for the corresponding small-polaron model 



have been constructed in the past 2CJ, |21[ . Studies of the spectral problem for these chains, 



however, have been restricted to the diagonal case. 
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APPENDIX A: SUPERALGEBRAS AND -MATRICES 



1. General linear Lie Superalgebras 

Let N, m,n G N and {e 3 - )c f Y a gl'i m+n be a homogeneous basis of an associative super- 
algebra, subject to the commutation relations 

bJWL = S jh (^-(-lyMMsie^) , (Al) 



whereas [X, Y} ± = XY - (-1) P ^ P ^YX denotes the so-called super commutator and p(X) 
gives the parity of a homogeneous element X of the superalgebra, that is 

/•(.Y)H ■ (A2) 

1 if X is an element of the odd subspace . 

Considering the super commutator as a generalized Lie product, the generators ej,£ consti- 
tute the Lie superalgebra gl(m\n). We restrict ourselves to the special case m = n = 1. By 
identifying 

Cj = e jn 2 , c) = ej,/ , n j = c]cj = e i)2 2 and n j = cf) = \-n- = e j: ^ (A3) 

we find gt(l|l) to be the algebra T of operators and creating and annihilating spinless 
fermions on a one-dimensional lattice respectively, j being the site index. In this case the 
even subspace is spanned by and n- while cj and span the odd subspace. 

For a more detailed introduction to the construction of superalgebras on graded vector 



spaces, we refer to [22j, [15j and section 12.3 in 



1 . 



2. Grassmann algebras 

Grassmann numbers, being the elements of a Grassmann algebra, are one of the key 
ingredients in the formulation of non-diagonal boundary conditions for super spin chains. 
The J\f 6 N generators of a Grassmann algebra will be denoted by £2, ■ ■ ■ , £m and in 



accordance with [22[] we define a product between them such that for all j, k, I = 1, 2, . . . , Af 



1. the product is associative, 

(E j S k )S l = £ j (S k S l ), (A4) 
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2. any two generators mutually anticommute, 

£j£ k = -£ k £j, (A5) 

3. and each non-zero product 

£ h £ h ...£ jr ,l<r<M (A6) 

involving r generators is linearly independent of products involving less than r gener- 
ators. In particular, this means that Grassmann generators £j have no inverse. 

For consistency reasons it is customary to supplement the set of generators by an iden- 
tity 1 with the defining properties 1-1 = 1 and l£j = £jl = £j. Using multi-index 
notation, each product of generators can be written as £ ll = £j 1 £j 2 ■ ■ whereas 
/i = {jijjz, ■ ■ ■ , is an, without loss of generality, ascendingly ordered set of natural 
numbers 1 < j n < M . The identity may be incorporated by setting £0 = 1. Finally, 
this enables us to express every Grassmann number G as a linear combination of generator 
products £ ii with complex coefficients G^, 

G = (A7) 

Here the summation is to be carried out over all multi-indices fi. In the following text 
this complex Grassmann algebra with M generators will be labeled CG^f. We impose a 
convenient grading, setting 

p{S tl ) = \n\ mod 2. (A8) 

The complex conjugation of a Grassmann number G is given by the complex conjugation of 
the linear coefficients in (IA7I) . i.e. G* = (G >1 )*£ IM . Moreover we define the adjoint G^ of a 
Grassmann number G by 

G « = (-i)*K^)(G^)*£ M . (A9) 



3. Super matrices 

Just like the elements of the above super algebras, super matrices are graded objects. Here 
we will only make use of square even invertable super matrices M , having the partitioning 

(a b\ 

M=\ , (A10) 
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such that all entries of the submatrices A and D are even elements of a superalgebra, whereas 
all entries of the submatrices B and C are odd elements of the same superalgebra. We define 
convenient analogs to the usual matrix operations. The supertrace is given by 



str{M} = tr{ A} - tr { £> } . (All) 

In contrast to the ordinary matrix transposition, the super transposition ( ) st is not an 
involution. Therefore, we have an additional inverse super transposition ( ) lst , 

t ( AT c T \ . f (a t -c t \ , 

M st = , M lst = . A12 

I -B T D T ) \B T D T J 

If the underlying superalgebra is CGjv" there are two more important operations, namely 
the adjoint operation 

Mt s ( {C 'A , (A13) 

\( B >r (D>r) 

where A^ is defined by entrywise application of ( 1A9I) . and the multiplication of a super matix 
by a Grassmann number G of definite parity, 

,Gl d w \ (a b\ 

G-M= . (A14) 

(-l) p(G) Gl dimD / \C D 
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